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The prospects for the measurement of 6^ oscillations with the ATLAS detector at the Large Hadron Collider are presented. B° candi- 
dates in the D^n^ and Z)7a| decay modes from semileptonic events were fully simulated and reconstructed, using a detailed detector 
description. The sensitivity and the expected accuracy for the measurement of the oscillation frequency were derived from unbinned 
maximum likelihood amplitude fits as functions of the integrated luminosity. A detailed treatment of the systematic uncertainties was 
performed. The dependence of the measurement sensitivity on various parameters was also evaluated. 



1 Introduction 

The observed and states are linear combinations of 
two mass eigenstates, denoted here as H and L. Due to the 
non-conservation of flavour in charged weak-current inter- 
actions, transitions between and B° states occur with a 
frequency proportional to Anig - mn - mi. 

Experimentally, the B'' oscillations have not yet been 
observed dkectly. The combined lower limit from mea- 
surements done by the ALEPH, DELPHI and OPAL exper- 
iments at LEP, by SLD at SLC, and by CDF at the Tevatron, 
is Am.5 > 14.4 ps"' at 95% CL, with a sensitivity at 95% 
CL of 19.3 ps ' [1]. In the Standard Model, it would be dif- 
ficult to accommodate values of Am^ above ~ 25 ps"' [2]. 

In this paper, the prospects of the ATLAS experiment at 
the Large Hadron Collider (LHC) to measure B^- B^l oscil- 
lations are presented. A detailed description of the analy- 
sis on which this presentation is based on can be found in 
Ref. [3]. A short discussion of subsequent changes is also 
included. 

2 Event selection 

The signal channels considered in this analysis for the 
measurement of B^- B'j oscillations are B'' — > DsTt and 
B^l — > Djfli, with Ds — > (pn followed by — > K^K , se- 
lected in semileptonic events. 

The event samples from this simulation study were gener- 
ated using PYTHIA 5.7 [4], passed then through a detailed 
GEANT3-based simulation of the ATLAS Inner Detector; 
charged tracks were then reconstructed using an algorithm 
based on the Kalman filter. The production of the bb- 
quark pairs in collisions at a centre-of-mass energy of 
= 14 TeV included direct production, gluon splitting, 
and flavour excitation processes. The fe-quark was forced 

'Speaker at the Workshop. 



to decay semileptonically giving a muon with transverse 
momentum pj > 6 GeV and pseudo-rapidity \ri\ < 2.5 
which is used by the level- 1 trigger to select the B hadronic 
channels in ATLAS, while the associated b was forced to 
produce the required B-decay channels. 

A multi-level trigger is used in ATLAS to select events. For 
B-physics, the level- 1 trigger is an inclusive muon trigger, 
as mentioned before. The level-2 trigger reconfirms the 
muon from level- 1 trigger using also the precision muon 
chambers, then in an un-guided search for tracks in the 
Inner Detector reconstructs a meson and, adding a new 
track, a meson. The level-3 trigger (the event filter) con- 
firms the level-2 result using a set of loose offline cuts to 
select the events. 

The flavour of the B^ meson, i.e. the particle or antiparticle 
state, is tagged at the production point by the muon used for 
the level- 1 trigger; at the decay vertex, the meson's state is 
given by the charge of the reconstructed Dj meson. 

Offline, the B*' meson was reconstructed from its de- 
cay products, applying kinematical cuts on reconstructed 
tracks, mass and vertex-fit cuts on the intermediate parti- 
cles, and cuts on properties of the B*' candidates (vertex-fit 
quality, proper time, impact parameter, mass, etc.). 

The background was estimated considering various four- or 
six-body B-hadron decay channels, and the combinatorial 
background. The four- and six-body background events 
were generated, passed through the detailed detector simu- 
lation program, reconstructed and analyzed using the same 
programs, the same conditions and the same cuts as the sig- 
nal events. For the study of the combinatorial background, 
about 1 . 1 million bb — > fiX events were analyzed using a 
fast detector simulation package. 

The expected number of signal and background events for 
an integrated luminosity of 10 fb"', corresponding to a 
one-year run at 10^^ cm^^s ' (so-called "low luminosity") 
are summarized in Table 1. To compute them, branching 
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ratios from Ref. [5] were taken, where known, else from 
PYTHIA. Charge-conjugate channels were taken into ac- 
count; corrections for trigger efficiencies of 63% (level-2 
Ds) and 82% (muons, level- 1 and offline combined) were 
applied. 





Process 


Events 


Signal 




2370 


channels 




870 






<400 


Exclusive 




<340 


background 


B° ^ D tt^ 


3 


channels 


B° D-a\ 


1 






2 









Combin. 


4 charged tracks 


1920 


background 


6 charged tracks 


1830 



Table 1. Signal and background samples analyzed for the study 
of fij- 6"J oscillations, and numbers of events expected for 10 
fb-'. 



The reconstructed 5° invariant-mass distribution in the de- 
cay channel ~* ^7^^ is shown in Figure 1 for an inte- 
grated luminosity of 10 fb~'. 



3 Extraction of significance limits and accu- 
racy for the measurement of Affij 

3.1 Proper-time reconstruction and resolution 

The proper time of the reconstructed B^ candidates was 
computed from the reconstructed transverse decay length, 
dxy, and from the Z?" transverse momentum, pj: 



t = = d_,yg 

cpT 



where g = Mgo/{cpj) and M^o is the Z?" mass. Its resolu- 
tion function Res(f | fo) was parameterized with the sum of 
two Gaussian functions, see Eq. (1), with parameters given 
in Table 2 for the signal channels. Here fo denotes the true 
(generated) proper time; /„ is the fraction and ctq, the width 
of the Gaussian function a. Similar parameterizations were 
obtained for the background channels Z?o — > D^n^ and 
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Figure 1. Reconstructed B\ invariant-mass distribution for B° 
D^^n* decays. Open histogram: signal, hatched: background from 
SjJ — » Dj^r^ decays, dark: fake reconstructed decays from the 
signal sample. Combinatorial background not shown here. 











D-a\ 


a 


fa (%) 


o-„(fs) 


fa (%) 


0-u,(fs) 


1 


59.6 + 6.6 


51.5 ±4.0 


62.5 + 14.1 


51.6+ 6.4 


2 


40.4 + 6.6 


107.3 ± 8.5 


37.5 + 14.1 


92.8 ± 12.7 



Table 2. Proper-time resolution function Res(f 1 1^) parameteriza- 
tion with the sum of two Gaussian functions. 



3.2 Likelihood function 



The probability density to observe an initial B*! meson {j - 

d, s) decaying at time to after its creation as a fiO meson is 
given by: 



Pj(to, yuo) 



2r, 



AE/o 

cosh h yU() COS AnijlQ 



(2) 



where AE^ 



E^ - E[, E,. = (E^ + E[)/2 and = -1. 
Eor the unmixed case (an initial 5° meson decays as a B^. 
meson at time Iq), the probabiUty density is given by the 
above expression with yUo - +1. 
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The above probability is modified by experimental ef- 
fects: finite proper-time resolution, wrong tags at produc- 
tion or decay, and background. Convolving p/fo, //q) 
with the proper time resolution ReSj(f | fo), one obtains 
the probability as a function of yUo and the reconstructed 
proper time t: qjit,fio) = J^^ Pj(fo,f^o) Res/f | fo) dfo, 
with a normalization factor and = 0.4 ps the cut 
on the proper decay time. Assuming a fraction ajj 
of wrong tags at production or decay, the probability 
becomes q'jit,fj.) = (1 - (i)j)qj{t,^) + Ujqj{t, -/u). Includ- 
ing the background, composed of oscillating B° mesons 
and of combinatorial background, with fractions fj (j - 
s, d, and combinatorial background cb), one obtains: 
pdfi(f,;u) = Y,j=s.d.chfj [(1 - u)i)qi{t,n) + (^jqjit, 
where the index ^ = 1 denotes the — > D^n^ chan- 
nel and = 2 the B',' 
the total sample is written as 



channel. The likelihood of 



N„ NL 



i:(Am,,Ar,) = ]~[ ]~[pdf,(f„;/,) 



(3) 



where A^^^ is the total number of events of type k, and 
A^ch = 2. 

3.3 Significance limits for the measurement of AaMj 

The ATLAS sensitivity for the Am^ measurement was de- 
termined using a simplified Monte-Carlo model to produce 
event samples, combined with the amplitude-fit method [6] 
to extract the limits. In the amplitude-fit method a new pa- 
rameter, the 5° oscillation amplitude y[, is introduced in 
the likelihood function by replacing the term '/io cos Aw sfo' 
with '/io-^cos Amyfo' in the B*' probability density func- 
tion. For each value of Am,, the new likelihood function 
is minimized with respect to yi, keeping all other parame- 
ters fixed, and a value ± o-*^"' is obtained. The statistical 
significance S of an oscillation signal can be expressed as 
S ~ 1 /cj^. One defines a 5cr significance limit as the value 
of Am.s for which l/crji - 5, and a sensitivity at 95% con- 
fidence limit as the value of Am.s for which I/ctj^ - 1.645. 

For Am, values smaller than the 5cr significance limit, the 
expected accuracy is estimated using the log-likelihood 
method, with the likelihood function given by Eq. (3). 



was defined, reducing the /^o error and the uncertainties on 
the widths from the proper time parameterization to values 
expected at the time of ATLAS data taking. 

4 Results and conclusions 

Table 3 shows the dependence of the amplitude and its 
statistical and systematic uncertainties on Am, for an inte- 
grated luminosity of 10 fb"', for both actual and projected 
systematic uncertainties. In the generated event samples, 
the value of Am, was set to Amf^" = oo, therefore one ex- 
pects cr^ compatible with zero. The dominant contribu- 
tions to the systematic uncertainty come from the uncer- 
tainty on the fgo fraction and from the parameterization of 
the proper time resolution. 



Am, ps ' 


10 ps-' 


20 ps-' 


30 ps-' 


y[ 0.045 


0.189 


0.042 


-0.291 


o-f +0.048 


+0.090 


±0.167 


±0.357 


syst +0.097 
" J\ -0.084 


+0.130 
-0.096 


+0.180 
-0.142 


+0.298 
-0.226 


with 'projected systematic uncertainties' 


syst +0.049 
" -0.049 


+0.060 
-0.048 


+0.085 
-0.066 


+0.137 
-0.117 



Table 3. The oscillation amplitude ^ and its statistical and sys- 
tematic uncertainties as a function of Am^ for an integrated lumi- 
nosity of 10 fb-'. 




3.4 Systematic uncertainties 

An attempt to estimate the systematic uncertainties was 
done. The following contributions to the systematic un- 
certainties were considered: a relative eiTor of 5% on the 
wrong-tag fraction for both B° and B^; ±lcr variation of 
Gaussian-function widths from Res(f | to) parameterization; 
fgo - BR(b — » B°), B° lifetime, Am^ varied separately 
by PDG uncertainty; 5% uncertainty for decay time r^b of 
combinatorial background, keeping the shape exponential. 
An additional set of 'projected systematic uncertainties' 



Figure 2. The significance of the S° oscillation signal as a func- 
tion of Am, for an integrated luminosity of 10 flj-' . 

The significance of the B^ oscillation signal as a function 
of A77ij for an integrated luminosity of 10 fb-'is shown in 
Fig. 2. The 5<j significance limit is 22.5 ps ' and the 95% 
CL sensitivity is 36.0 ps-', when computed with the statis- 
tical uncertainty only. Computed with the total uncertainty, 
the 5cr significance limit is 16.0 ps ' and the 95% CL sen- 
sitivity is 34.5 ps-' for the actual systematic uncertainties. 
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and 21 ps ' and 35.5 ps ' for the projected systematic un- 
certainties. The hmits for other values of the integrated lu- 
minosity are given in Table 4, computed with the statistical 
uncertainties only. 



Luminosity 


5cr limit 


95% CL sensitivity 


(ft-') 




(ps-') 


5 


17.5 


32.0 


10 


22.5 


36.0 


20 


27.0 


39.0 


30 


29.5 


41.0 



Table 4. The dependence of significance limits for the Amj mea- 
surement on the integrated luminosity, computed with statistical 
uncertainties only. 



The dependence of the significance limits on ATj/Ts was 
also estimated. For values of ATj/Fj up to 30%, no size- 
able effect was observed. The shape and the fraction of 
the combinatorial background were also varied within rea- 
sonable values; only a weak dependence of the limits was 
observed. 

For Am, values smaller than the 5<t significance limit, the 
accuracy of the Am, measurement was determined for dif- 
ferent values of the integrated luminosity. The results are 
given in Table 5. If measured, the precision on Am, will be 
dominated by the statistical errors. 



Luminosity 


Amf™ 


A rec , Am. , Am. 


Obs. 


(ft-') 


(ps-') 


(ps-') 




5 


17.5 


17.689 ±0.083 ±0.002 


5crlim. 


10 


15.0 


L5.021 ±0.049 ±0.002 




22.5 


22.396 ±0.072 ±0.005 


5(r lim. 




15.0 


14.949 ± 0.033 ± 0.002 




20 


20.0 


20.041 ±0.068 ±0.005 






27.0 


26.948 ± 0.070 ± 0.003 


5cr lim. 




15.0 


14.942 ± 0.028 ± 0.004 




30 


20.0 


20.010 + 0.043 + 0.002 






29.5 


29.708 ± 0.083 ± 0.007 


5crlim. 



Table 5. The accuracy of Awj^ measurement as a function of the 

integrated luminosity, cr^^,' represents the statistical uncertainty, 
crf^' the systematic uncertainty. The 5cr limits were computed 
with statistical errors only. 



5 Recent developments 

This section summarizes recent changes of the assumptions 
or conditions used to get the results presented above. 

The most important changes in the detector geometry are 
the increase of the beam-pipe diameter from 41.5 mm to 
50.5 mm and the increase of the pixel length in B-layer, the 
closest layer to the beam-pipe, from 300 //m to 400 jum. 



Due to financial constraints, the B-physics trigger re- 
sources have to be minimized. Previously, dedicated re- 
sources were supposed to be available for the B-physics 
trigger, in addition to resources for high-/?T physics trigger 
('discovery physics trigger'). It may not be possible to pro- 
vide any significant additional resources. Moreover, there 
are financial uncertainties which could lead to the deferral 
of some detector items, therefore it is possible to have a 
reduced detector at start-up. Items included in the deferral 
scenario are the second pixel layer from the Inner Detector, 
resulting in a two-layer pixel detector, the 77 > 2 region of 
the Transition Radiation Detector from the Inner Detector, 
and a significant part of the processors for level-2 and event 
filter, reducing the computing resources for the high-level 
trigger and limiting the level- 1 rate. 

The luminosity target for LHC start-up doubled to 2 x 
10^^ cm"^s"\ therefore it will be necessary to re-evaluate 
trigger thresholds and to remove some triggers requiring 
too much resources. The trigger for oscillation chan- 
nels requires a significant rate, therefore it is very likely 
that the muon trigger threshold will be raised to ~ 8 GeV. 

Recent work concentrates on trigger-related issues; im- 
provements of the offline analysis, although possible, have 
to be postponed. To reduce the resource requirements, one 
of the possibilities would be to change the trigger for B 
hadronic channels from juipj > 6 GeV) at level- 1 plus 
Inner Detector full scan at Level-2 to luip-r > 6 GeV) 
at level- 1, plus a low Et level- 1 calorimeter Region-of- 
Interest (Rol), used then to guide reconstruction at level-2. 
In addition, one can use the level-2 Rol to limit the region 
for reconstruction at the event filter. A flexible trigger and 
analysis strategy is nowadays evaluated, which should be 
able to cope with detector changes, luminosity scenarios, 
financial problems. Another direction is to recover perfor- 
mance using optimized reconstruction algorithms, flexible 
trigger and analysis thresholds. A re-evaluation of the de- 
tector performance with the latest geometry is also under 
way. Preliminary results are encouraging. 
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